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optimal control theory, apply seta rule for spacecraft transcribe the trajectory 
Euler-Lagrange theorem and _ steering — a suboptimal optimization into a 
solve two point boundary policy that is acceptable by parameter optimization 
value problem (TPBVP) the mission designer problem 
Edelbaum Kluever Betts 
Alfano Chang Whiffen 
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1 enh one HDDP 
¢ Hybrid Differential Dynamic Programming 
¢ introduced by Lantoine and Russell 


* sequence of control updates that minimize quadratic model of cost-to-go 
* map derivatives along trajectory with state transition matrix and tensor 


19 

\ 
Ni 
‘e | 


1/4/2017 


Celestial and Spaceflight | DDP 
Mechanics Laboratory 


¢ Hybrid Differential Dynamic Programming 
¢ introduced by Lantoine and Russell 
* sequence of control updates that minimize quadratic model of cost-to-go 
* map derivatives along trajectory with state transition matrix and tensor 


Quadratic model 
from derivatives of w 


1/4/2017 


Celestial and Spaceflight | DDP 
Mechanics Laboratory 


¢ Hybrid Differential Dynamic Programming 
¢ introduced by Lantoine and Russell 
* sequence of control updates that minimize quadratic model of cost-to-go 
* map derivatives along trajectory with state transition matrix and tensor 


Quadratic model 
from derivatives of w 


v 
ata se. 
at 
. ~~. 
'. s, 
%e 
*, 


Derivatives mapped 
backward with STMs 


1/4/2017 


Celestial and Spaceflight | DDP 
Mechanics Laboratory 


¢ Hybrid Differential Dynamic Programming 
¢ introduced by Lantoine and Russell 
* sequence of control updates that minimize quadratic model of cost-to-go 
* map derivatives along trajectory with state transition matrix and tensor 


Quadratic model 
Minimize local model from derivatives of 


Uy PT ae U3 OU3 a “ . 


v 
ata se. 
ay te 
« ~s, 
e 
%e 
*, 


Derivatives mapped 
backward with STMs 


1/4/2017 


Celestial and Spaceflight | DDP 
Mechanics Laboratory 


¢ Hybrid Differential Dynamic Programming 
¢ introduced by Lantoine and Russell 
* sequence of control updates that minimize quadratic model of cost-to-go 
* map derivatives along trajectory with state transition matrix and tensor 


Quadratic model 


Minimize local model from derivatives of 


“Tre 
To COT nee 


ane 
a, ns 


°. 
%e 


oo are 
tes Ve 
poe 
oe 
2 


Derivatives mapped 
backward with STMs 


1/4/2017 


Celestial and Spaceflight | DDP 
Mechanics Laboratory 
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¢ Forward pass: evaluate u + du in equations of motion 
¢ Backward sweep: compute each du; 
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See: 

¢ Forward pass: evaluate u + du i Gregory Lantoine and Ryan P. Russell. A hybrid differential 
dynamic programming algorithm 

° Backward sweep: compute each for constrained optimal control problems. part 1: Theory. 
Journal of Optimization Theory 

and Applications, 154(2):382-417, 2012. 


1/4/2017 au 


oy 


aries pul prifeant el The Sundman Transformation 


¢ Change independent variable from time to a function of orbital radius 
dt = cyr"dt 


¢ Can choose n, Cy, so that T is an orbit angle 


Eccentric Anomaly Mean Anomaly True Anomaly 


a a? r? 
dt = |-rdE Apo dt = —dv 
be Le h 


¢ Equations of motion become x’ = xcy,r” ¢ Track time in the state vector 
¢ Discretize in T : : 
x x 
¢ Specify To, Ty, rather than to, ty x=|y], x’ = lylcyr” 
Z 


* 1.e. specify number of revolutions 
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HDDP + Sundman 


¢ Choose T = E, the eccentric anomaly 
¢ Minimum fuel GTO to GEO in 450.5 revs 
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1 rename HDDP + Sundman 


¢ Choose T = E, the eccentric anomaly 

¢ Minimum fuel GTO to GEO in 450.5 revs 
¢ 2-body dynamics + J, 

135,150 variables 


¢ 61 minutes 
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HDDP + Sundman 


¢ Choose T = E, the eccentric anomaly 

¢ Minimum fuel GTO to GEO in 450.5 revs 
¢ 2-body dynamics + J2+ lunar gravity 
135,150 variables 


¢ 107 minutes 
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1 rename HDDP + Sundman 


¢ Choose T = E, the eccentric anomaly 

¢ Minimum fuel GTO to GEO in 1000.5 revs 
¢ 2-body dynamics + J2+ lunar gravity 
300,150 variables 
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Table 2.1: Summary of GTO to GEO Results. 
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500 rev orbit lowering at Mars 
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1 reno Trust-Region Quadratic Subproblem 


¢ Feedback control laws for du, 6A, dw are unconstrained 

¢ Likely to step beyond validity of quadratic expansion 

¢ Require invertible, positive definite Hessians (negative definite for J), ) 
¢ Trust-region quadratic subproblem (TRQP): 


L. : 
min| J; dup + ~dUj Juu~OUr| 
OU i 2 
s.t. || Dduz|| < A 
¢ Acceptance of an iterate: 


5] 7 min((1+#«)A»y,Amax), if o€ [1—e1,14+ 41] 
= p+1— 
"~ ERoo - 


max((1—#)A,, Amin), otherwise 
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